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NOLIN 

A NONLINEAR LAMINATE ANALYSIS PROGRAM 

by 

John J. Kibler 


SUMMARY 

A nonlinear, plane-stress, laminate -analysis program, 
NOLIN, has been developed which accounts for laminae non- 
linearity under inplane shear and transverse extensional 
stress. The program determines the nonlinear stress-strain 
behavior of symmetric laminates subjected to any combination 
of inplane shear and biaxial extensional loadings. The pro- 
gram has the ability to treat different stress-strain be- 
havior in tension and compression, and predicts laminate 
failure using any or all of maximum stress, maximum strain, 
and quadratic interaction failure criteria. 

A second program, UNI, has been developed which com- 
putes elastic constants and thermal coefficients of expan- 
sion for laminae, from constituent properties, to aid in 
compiling input for the NOLIN program. Laminae properties 
can be computed for isotropic or transversely isotropic 
fibers in an isotropic matrix. In addition, the nonlinear 
inplane shear stress-strain curves are computed for the 
laminae by computing the Ramberg-Osgood shear stress para- 
meter. 

This document provides brief descriptions of both 
programs , a description of the flow of information through 
the NOLIN program, and detailed descriptions of the input 
required for each program. Sections are provided with 
sample problems and sample program output , along with 
complete listings of each program. 



1. NOLIN Program Description 
1.1 Introduction 

The NOLIN program codifies a nonlinear, plane-stress, 
laminate analysis wherein the nonlinear behavior of the 
laminae under inplane shear and transverse extensional 
stresses are taken into account. Both the underlying analy- 
tical development and the computer program are sufficiently 
general to enable the user to study the nonlinear behavior 
of a symmetric laminate subjected to any combination of in- 
plane shear and biaxial extensional, loadings. Contained in 
this document are a detailed description of the input re- 
quired to use the NOLIN program, as well as a description of 
the theoretical developments upon which the program is based. 
For a complete theoretical description the reader is referred 
to Re f . 1 . 

In unidirectional, fiber-reinforced laminae, the trans- 
verse extensional and, particularly, the inplane shear stress- 
strain relationships cannot be accurately characterized as 
linear. The NOLIN program allows for nonlinear representations 
by permitting these stress-strain relationships to take the 
form of Ramberg-Osgood nonlinear relationships. The introduc- 
tion of these nonlinear, Ramberg-Osgood type constitutive 
relationships into a laminate analysis then leads to a set of 
nonlinear equations involving the laminae stress components 
as unknowns . The program then solves this set of equations 
by means of a generalized, Newton-Raphson procedure to give 
the laminae stresses and strains corresponding to the applied 
boundary stresses. 

The theoretical development for this nonlinear, laminate 
analysis incorporates total deformation theory with Ramberg- 
Osgood type stress-strain characterizations to formulate the 
governing nonlinear equations. At the outset the compliance 
tensor is assumed to be the sum of two tensors, the components 
of one are the usual components associated with linear, 
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orthotropic, plane-stress elasticity theory, while the 
second tensor contains the nonlinear elements. By assuming 
a quadratic interaction of the stress components, and re- 
quiring the constitutive relationship to reduce to the re- 
lationships for the uniaxial stress cases of inplane shear 
and transverse extension, the elements of the nonlinear com- 
pliance tensor are explicitly determined. 

Having the nonlinear laminae constitutive relations , 
the usual methods of laminate theory are then utilized to 
obtain the governing, nonlinear equations for the laminate. 

As in linear, laminate theory, the strains of the individual 
laminae are first rotated to a common set of laminate axis, 
and the laminate compatibility relations requiring tl\e corres- 
ponding strains of the individual laminae to be equal are 
then employed. In addition, equilibrium at the laminate 
boundaries is envoked. In this way the required number of 
equations involving the unknown laminae stresses are formulated. 

The program solution procedure for the set. of nonlinear 
equations involving the laminae stress components is a Newton- 
Raphson technique generalized to accomodate systems of equa- 
tions. The starting point for the solution procedure is taken 
as the solution of the associated, linear laminate problem, 
where the associated linear problem is obtained by ignoring 
all nonlinear terms. 

Incorporated into the program are three different failure 
criteria, maximum stress, maximum strain and a quadratic inter- 
action criteria. Any or all of these may be employed by the 
program user. Unfortunately, the nonlinear aspects of this 
program preclude the generation of strength envelopes since 
linear extrapolation is not valid here. Instead a sequence 
of combined loadings may be run. 
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1.2 Derivation of Governing Equations 

In the usual , plane-stress , linear, laminate analysis, 
the laminae stress-strain relationships are first written 
with respect to the longitudinal and transverse fiber direc- 
tions. These laminae relationships are then transformed to 
a common set of laminate axes and the laminate constitutive 
relationships are then developed by enforcing strain compa- 
tibility. In this nonlinear laminate analysis an entirely 
analogous procedure is followed. The laminae nonlinear, 
stress-strain relationships are first written for the N 
lamina with reference to their respective natural axes, and 
then transformed to a common set of laminate axes. Enforce- 
ment of compatibility between adjacent laminae as well as 
equilibrium at the laminate boundaries leads to the govern- 
ing, nonlinear equations. 

Let the strains be separable into linear and nonlinear 
parts , 
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Where primes designate the linear parts and the double primes 
designate the nonlinear parts, and the subscripts L and T 
refer to the longitudinal and transverse directions, respect- 
ively. The linear parts of the strain components are related 
to the stress components by the usual orthotropic, plane- 
stress elasticity constitutive relationship, 
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For the uniaxial stress states of transverse extension 
and inplane shear,, the nonlinear parts of the strain com- 
ponents are related to the stress components by the re- 
lationships, , 
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Where the values of n, m, f^ and are obtained from 
curve fitting unidirectional, experimental data. It is 
assumed that in the longitudinal direction, the nonlinear 
part of the strain component is zero. 

i i 


The relationships between the nonlinear parts of the strain 
components and the stress components for cases of combined 
transverse-extension and inplane-shear stresses are deter- 
mined by using a generalization of ^ theory in addition to 
the requirement that these relationships teduce to equations 
(3) for uniaxial states of stress. The resulting equations 
are , 


n-1 



Equations (2) and (5) comprise the total, nonlinear, stress- 
strain relationship for the laminae. 


5 



The governing equations are formulated so that the 
three stress components in each lamina are the unknowns. Thus, 
for an N-layered laminate, the problem is formulated in terms 
of 3N unknowns. To obtain solutions, 3N equations are then 
required, and these equations consist of three equilibrium 
equations and 3(N-1) compatibility equations satisfying strain 
compatibility between adjacent laminae. The three equations 
of equilibrium for a laminate under a combined state of stress 
are, 
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Where N.^, ^2 and are t ' ie a PPli et ^ stress resultants, t^ 

the thickness of the kth lamina, and subscripts 1 and 2 de- 
note the laminate axes. The 3(N-1) equations of strain 
compatibility are. 
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k = 2 , 3 , . . . , N 

Equations (6) and (7) are the 3N equations required for the 
solution of the nonlinear laminate problem. When' the stress- 
strain relations given by equations (2) , (4) and (5) are 

transformed to the laminate reference axes and substituted 
into equations (7) , the governing equations can be expressed 
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in functional form as , 


( o^ , o 2 > • • • / o ^ , •••) — 0 

k = 1, 2, . . . , 3N 
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1.3 Method of Solution 


Solutions of equations (8) for the 3N stress components 
are obtained by employing a Newton- Raphson iterative scheme. 
The functions F. are first expanded in Taylor series about 

K 1 0 

an approximate set of initial stresses, a ^ . Considering 
only the first order terms of these series, 
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For clarity, the notation in equations (10) is, in expanded 
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The solution for in equation (10) may be taken as the 
approximate, initial stress values for the next iteration 
step, and this process repeated until a result is obtained 
within some desired accuracy. After the stresses are ob- 
tained and transformed to the laminae natural axes, the 
corresponding laminae strains are determined from equations 
(1) , (2) and (5) . 
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1 . 4 Computer Program 

The flow chart for the computer program is shown in 
Fig. l. The major sections of the program are the forma- 
tion of the governing equations, the Newton- Raphson solu- 
tion procedure and the failure checks. 

Using the computer program. notation , the governing 
equations take the form, 

[A] • SG + B = SGO (15) 

Where SG and SGO are the stress solution vector and the 
applied stress vector, respectively. A is a matrix of con- 
stant elements which are the coefficients of the linear 
terms in the solution, and B is a vector containing the 
nonlinear terms in the solution. The set of equations (15) 
are equivalent to equations (8) . If in equation (15) the 
vector B is set to zero, the resulting equation 

[A] * SG = SGO (16) 

is the linear laminate solution. The stress vector, SG, 
as determined from equation (16) is taken as the initial 
approximation for the stress vector in the Newton-Raphson 
procedure. 

For the Newton-Raphson procedure it is necessary to 
formulate the derivative of ([A]* SG + B) in equation (15) 
with respect to a ^ as well as the vector 

DC = ( [A] *SGO + B - SGO), (17) 

The vector DC corresponds to the vector Fjf in equation 

(10), and an explicit evaluation of DC is obtained by 

using the current, approximate value for the solution 

stress vector, SC. The derivative of ([A]*SG + B) is 

designated DB in the computer program, and is equivalent 

to the matrix (9F,/3a.) a? in equation (10). An explicit 
' K 1 1 

evaluation of DB is also, obtained by using the current, 
a PProximate value for the solution stress vector, SG. 
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In the program, the external loading is applied in. incre- 
ments. The approximate solution stress vector for the first 
load increment and the first Newton- Raphson iteration is deter- 
mined from equation (16) . For the second and third load in- 
crements, the approximate solution stress vectors for the first 
iteration are taken as the final solution stress vectors from 
the previous increments. Solutions for subsequent load in- 
crements are initiated by the following algorithm: 


< SG i + 1 ^ INITIAL - (SG i>FIHAL * <FACT0E > 

(FACTOR) - LJiZ$. l^i’riBAL ~ (SG i-l>FraM, 

(i-ir b l 'FINAL 


(18) 


The convergence and divergence criteria employed in the 
program are contained in the following expressions: 

| (SG i+1 - SG i )/SG i | < e 

(19) 

| (SG i+ i - SGp/SGj < X 


Where SG. and SG.,. are the solution vectors obtained from 

the i and i+1 1 " iterations. Usually values of 

- 3 +4 

10 and 10 are taken for e and A, respectively. ' However, 

• V. * 

the other values may be input as data to the program. In 
addition, the maximum number of iterations to be allowed is 
input as data. Ten iterations have been found to be suffi- 
cient for most problems. 

The program contains three failure criteria, maximum 
strain, maximum stress, and a quadratic interaction criteria. 
After a solution is obtained for each load increment any or 
all of these failure criterias may be applied to check for 
laminae failure. 

The maximum stress and maximum strain failure criteria 
check, respectively, the laminae stress or strain values in 
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the fiber, and transverse fiber directions against the 
material allowables. These allowables are input to the 
program as data. The quadratic criteria is given by 


A 11 °LL + A 22 °TT + A 44 a LT 

( 20 ) 

+A 12 °LL °TT + B 1 °LL + B 2 a TT = 1 

where the coefficients are functions of the allowable 
stress 



t C ' . 

F t and F t are the allowable tension and compression 
L L t c 

stresses in the longitudinal direction, F T and F T are the 

allowable tension and compression stresses in the transverse 
direction, and F s is the allowable shear stress. The co- 
efficient A ^2 is input as data to the program, or a default 
null value is used. 

If a failure criteria is satisfied at the end of a load 
increment, the program determines the failure load through 
linear interpolation. If all failure criteria are being 
checked, and not all indicate failure during the same load 
increment, the program continues loading until all criteria 
indicate failure. 
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Figure 1 


Solution flow Chart 














2. NOLIN PROGRAM USERS GUIDE 


2.1 Program Description 

This section describes the input data requirements 
for NOLIN (version 2 mod 2) . Input procedures have been 
streamlined wherever possible by using one NAMELIST state- 
ment, hence eliminating any input under format control. 

The input flow diagram and input description provide 
all information necessary to specify input data sets cap- 
able of exercising all program options. 

2.2 Input Description 

The initial data required is a message of five cards 
of alphanumeric descriptive information describing the 
problem being solved and printed as a title on the output. 
These five cards may be left blank, but must be included 
ahead of the first NAMELIST deck in the data. This descrip- 
tive message is read only once at the beginning of the pro- 
gram exectuion. The multiple case feature of running 
successive computations is accomplished by supplying multiple 
NAMELIST data sets with the changed variables indicated. 

The following is a description of the input variables 
required for execution of the program. Where appropriate, 
default or suggested values are indicated. The following 
data are supplied through NAMELIST "DATA" : 

Program Option Parameters 

IOPT: Ramberg-Osgood parameter sentinel. 

I0PT=1 : Input R-0 parameter directly. 
I0PT=2 : Determine parameters from stress- 
strain curve-fit routine. 

COPT: Curve-fit sentinel, (exercised if I0PT=2) 

C0PT=1 : Same stress-strain data for each 
layer. 

C0PT=2 : Different sets of data for each 


layer . 
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EOPT: 


Exponent option, (exercised if I0PT=2) 
E0PT=1: Determine exponents from curve- 
fit routine. 

E0PT=2 : Input exponents to curve-fit 
routine. 


Solution Accuracy Parameters 


KSGM: 
SMLT : 


EPS : 
UPBD : 
INMT: 


No. of load incrememts, maximum = 50 
increments. 

Load increment multiple , maximum number 
of Newton-Raphson iterations, default 
is 100 . 

Convergence criteria for Newton-Raphson 

-3 

analysis, default value is 10 
Divergence criteria during Newton-Raph- 
son analysis, default value is 20000. 
Incrementation estimate method, default 
value is 2. 


Layer Description 


NLAY: 

THICK (LAY) : 
THETA (LAY) : 
MATYPE(LAY) : 


Number of laminate layers (max. is 20). 
Thickness of each layer. 

Orientation of each layer in degrees. 
Material kind of each layer (maximum 
number of different materials is 20) . 


Material Description 


Ell (MATYPE) : 

E22 

G12 

VI 2 

SllT 

sue 

S22T 


Lamina longitudinal modulus. 

Lamina transverse modulus . 

Lamina shear modulus. 

Lamina major Poisson’s ratio. 

Lamina longitudinal tensile strength. 
Lamina longitudinal compressive strength 
Lamina transverse tensile strength. 
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S22C (MATYPE) : 

S12 

EP11T 

EP11C 

EP22T " 

EP22C 
GAMA " 

A12 


STY 

SCY 

TY 

XM 

XN 


Lamina transverse compressive strength. 
Lamina in-plane shear strength. 

Lamina longitudinal tensile strain. 
Lamina longitudinal compressive strain. 
Lamina transverse tensile strain. 

Lamina transverse compressive strain. 
Lamina in-plane shear strain. 

Lamina interaction term for quadratic 
interaction criteria - default value 
is 0.0. 

Ramberg-Osgood tension constant 
Ramberg-Osgood compression constant 
Ramberg-Osgood shear constant. 
Ramberg-Osgood shear exponent 
(default value is 3.0) . 

Ramberg-Osgood tension exponent 
(default value is 3.0). 


Stress-Strain Data (input in IOPT=2) 


IPTS : 

SIGH (I, MATYPE) : 

SIG22 

SIG12 

EPS11 

EPS22 " 

EPS12 


Number of stress-strain data points. 
IPTS values of longitudinal lamina 
stresses for each material type. 

IPTS values to transverse lamina 
stresses for each material type. 

IPTS values of in-plane shear stresses 
for each material type. 

IPTS values of longitudinal lamina 
strains for each material type. 

IPTS values of transverse lamina 
strains for each material type. 

IPTS values of in-plane shear strain 
for each material type. 
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Applied Loading and Failure Criteria 


SOU: 

S022: 

S012 : 
IFCN: 


STIFF : 


Initial axial stress applied to laminate 
Initial transverse stress applied to 
laminate . 

Initial shear stress applied to laminate 

Failure criteria sentinel . . 

IFCN=1: ultimate stress 

IFCN=2 : ultimate strain 

IFCN=3 : quadratic interaction 

IFCN=4 : all failure criteria 

Ratio of final to initial laminate 

stiffness which constitutes failure 

due to stiffness reduction, default 

value is 0.10. 
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3. UNI PROGRAM 


3.1 Description 

Program UNI computes the elastic properties , thermal 
expansion coefficients and the Ramberg-Osgood shear stress 
parameter for the unidirectional fiber bundle or the lamina. 
The fiber may be isotropic or transversely isotropic and the 
matrix is isotropic. The effective elastic properties of 
the composite are calculated from the composite cylinder 
assemblage model proposed by Hashin and Rosen [2] and the 
thermal expansion coefficients from the analytical results 
of Ref. [3] . The program also calculates the Ramberg-Osgood 
shear stress parameter of the matrix as outlined in Ref. [1] . 

The UNI program takes constituent properties as input 
and computes laminae properties as output. The output 
properties provide all the information required as input to 
the NOLIN program. The program accepts families of fiber 
and matrix materials such that an array of laminae proper- 
ties can be generated. This feature can be especially 
useful for sensitivity analyses. All input to UNI is 
accomplished through a single NAMELIST statement "UNID". 

3.2 Input Description 

The following describes the variables required for the 
execution of UNI through NAMELIST "UNID". 

Program Control Variables 

NF: No. of fibers, max. is 20. 

NM: No. of matrices, max. is 20. 

NVM: No. of matrix vol. fractions, max. is 20. 
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Matrix Properties, J=1 , NM 


EM (J) : 
RHOM(J) : 
ANUM(J) : 
ALPM(J) : 

ROMS ( J) : 


Young's modulus for Jth matrix. 
Density of Jth matrix. 

Poisson ratio for Jth matrix. 

Coef. of thermal expansion for Jth 
matrix required only if ISTS-1. 

Shear stress Ramberg-Osgood parameter 
for Jth matrix required only if 
NONLIN =1. . 


Isotropic Fiber Properties, J=l, NF 

The following variables are required when IS0T=1: 
EF(J): Young's modulus for Jth fiber. 


ANUF (J) 
RHOF (J) 
ALPF (J) 


Poisson ratio for Jth fiber. 
Density of Jth. -fiber. 

Coef. of thermal expansion for Jth 
fiber (required only if ISTS=1) . 


Transversely Isotropic Fiber Properties, J=1 NF 

The following variables are required when IS0T=2 : 


EFA(J) : . 
EFT (J) : 
ANUFA(J) : 
GFA(J) : 
ANUFT(J) : 
ALPF(J) : 

ALPFT(J) : 


Axial Young' s modulus for Jth fiber. 
Transverse Young's modulus for Jth fiber. 
Axial Poisson ratio for Jth fiber. 

Axial shear modulus for Jth fiber. 

Transverse Poisson ratio for Jth fiber. 

Axial thermal expan. coef. for Jth 
fiber (required only if ISTS=1) . . 

Transverse thermal exp. coef. for Jth fiber, 
(required only if ISTS=1) . 
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Laminae Volume Fractions, J=l, NVM 


VM(J) : 


Jth volume fraction of matrix material. 


The UNI program computes laminae properties for all 
combinations of fibers, matrix materials, and volume frac- 
tions which are supplied as input. That is, there are a 
total of NF x NM x NVM materials formed from the input. With- 
in a given run the total number of materials (NFxNMxNVM) 
must be less than 200. 

The following is a description of the properties com- 
puted by UNI and printed as output: 


Calculated Thermo-Elastic Constituent Parameters 


GFT ( J) : 
GF ( J) : 
GM'(J) : 
AKF(J) : 
AKM(J) : 


Transverse shear modulus for Jth fiber. 
Shear modulus for Jth fiber. 

Shear mdoulus for Jth matrix. 

Plane strain bulk modulus for Jth fiber. 
Plane strain bulk modulus for Jth matrix. 


Effective Thermo-Elastic Parameters 


AKTS (J) : 
EAS (J) : 
ETS (J) : 
ANUAS (J) : 
ANUTS (J) : 
GAS (J) : 
GTS ( J) : 


Effective trans. bulk modulus for Jth 
material. 

Effective axial Young's modulus for Jth 
material. 

Effective trans. Young's modulus for Jth 
material. 

Eff. Poisson ratio (unidirectional axial 
stress) for Jth material. 

Eff. Poisson ratio (in transverse plane) 
for Jth material. 

Eff. shear modulus (in fiber planes) for 
Jth material. 

Eff. shear mdoulus (in trans. planes) 
for Jth material. 
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ALP AS (J) : 


ALPTS { J) : 

RHOS (J) : . 
ROCOMP (J) : 


Eff. (fiber direction) thermal exp. 
coef. for Jth material. 

Eff. (trans. direction) thermal exp. coef. 
for Jth material. 

Bulk density for Jth material. 
Ramberg-Osgood shear stress parameter 
for Jth material. 
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4. SAMPLE PROBLEMS 


The purpose of this section is to present several sample 
problems which illustrate the capabilities of both the UNI 
and the NOLIN programs. Both the program input and output 
are listed at the end of this section to aid in understanding 
the program. 

4.1 UNI - Sample Problems 

4.1.1 Laminae properties for Thornel 50 fibers in two carbon 
matrices have been determined. The transversely isotropic 
fiber option has been used to model the Thornel 50 fibers, and 
a Ramberg-Osgood shear stress parameter is included for the 
two matrix materials. 

The input data and the resulting output data are given 
in sections 4.1.2 and 4.1.3. Note that the constituent proper- 
ties are mirrored in the program output along with the required 
computed constituent properties. The effective thermoelastic 
properties for the 1-D laminae are printed on the second page 
of output. Units are consistent throughout the program such 

that the units need only be consistent for the constituent 

2 

properties, in this case properties were input as MN/M . 

The second sample problem combines KEVLAR-49 fibers with 
a range of matrix properties to obtain the 1-D composite prop- 
erties. A list of input data for this case and the correspond- 
ing program output are shown in sections 4.1.2 and 4.1.3. It 
is interesting to note that an order of magnitude change in 
matrix modulus results in a factor of five change in trans- 
verse modulus and transverse shear modulus. 
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4.1.2 UNI SAMPLE PROBLEMS 


INPUT CARDS 


c un i n r a 


$ d a t aoup 

N F ~ 1 • ''i - - > , -a V ■' - . , 

EfA(i) S ?. ;•> t + o * ; • f > ( 1 ) = '/ . / "iK- + n «* , U F A 1 • = . 3 fi E ♦ (j 4 , 

A M I J F A ( 1 ) - t| . .t » A N f ! ( i. i s (i . 1. , 

R H 0 F (;)=.)., 

EM(1)=1. 7t:+. |4, • 

A N i.)M (1 ) = 2*.i . 2 * 

RHU U £ 1;)s'2*l , ; ; • • 

VM(l-)s-0.4,J}.6, 
rom's ( i ) =?*.s. n» 

A L P F ( 1 ,'F.-i )7, 

A L p F r .< 1 ) - 3 . o E - ij ■' , 

A L P M ( ) = 'l ■» 1. . 4 f: - i '■ > 
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ALP M ( I ) = !,!'*! . 39H- -r'5 • • 
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UNI SAMPLE PROBLEM OUTPUT 
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PROBLEM NO. 2 OUTPUT 
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EFFECTIVE thcRmO-ELASTIC P»R»MPTFflS 




rv « 

n >0 

n « 

rv « 

rv 4 

rv 0 

rv t* 

rv >o 

*> 0 

<* * 

O 4 

o. ca 

«p .«> 

* «j 

♦ 40 

<r « 


45 «■ 

« 0 

« 9 


• • 

• t 

• 1 

• # 

1 1 

• ■ 

• • 

• I 

1 • 

t • 

4 *- 

UJ ILt 

ai a> 

uj at 

uj uj 

tu Ul 

UJ ill 

UJ Ul 

Ul UJ 

UJ UJ 

UJ Uj 

w w 

a <M 

CV H 

i*v rv 

C3 >o 

0 . TS 

CM O 

-0 'O 

w in 

rv c\* 

0 0 

4 4 

O' CM 

'O <r 

90 4> 

0 m 

CM O 

in 0 

in 4 

4 ro 

c* 0 

in ro 

X X 

CO O 

o <r 

cm rv 

•0 a) 

to rv 

<-» M- 

H O 

n 4 

rv <© 

cm rv 

a. a 

ro to 

in cm 

a w 

n> c=» 

0 0 

C= <D 

tr N 

«d m 

cm ro 

rv w< 

.j _j 

rv ro 

PO 4- 

0 * tn 

0 <0 

m « 

0 rv 

in id 

T-| O 

cO 0 

tr n 

-« 4 


• m 

• V 

• V 

* — 

. • * 

• * 

• » 

— « 

• * 


ITk >0 

in <0 

4 <0 

V 0 

4 « 

»o <3 

fO « 

10 tO 

cm rv 

cm rv 


r 

t 

t . 

r 

• 

1 

1 

1 

1 

to 


« 

« * 

CD O 

4 O 

tO O 

rv 0 

to O 

in 0 

in 0 

*-♦ O 

<D O 

O O 

t- O 

- O 

• O 

* c 

• 0 

•O 

* 0 

- c 

• 0 

* O 

- O 

— X 

O O 

O 0 

O' O 

<n a 

tv O 

XI 0 

rv 0 

ro ci 

O O 

CM O 

x 0: 

O 

O C 

er 

CM O 

in 0 

fw CD 

4 c. 

ro cd 

ec cd 

CM O 


0 0 

rv 0 

v 0 

O a 

ro c? 

rv 0 

0 0 

ro 0 

in cd 

co a 


tr- • 

• 

«v • 

CM • 

ro • 

* ro v 

4 * 

•v * 

9 - 

4 • 


’ H 

H 


r* 


H 


H 

H 

rl 


♦ « 




















*-* — 

in 9 


ro 0 


V co 


CM 


«d vr 


M CC 


cm m 


CO t© 


'O 00 


rv ro 

-a >* 

to tC 


O - 


M - 


• • 


• m 


• * 


• • 


- # 


• • 


• • 

<^r 

1 «; • 


• - *-i 


• CM 


0 CM 


O CV» 


T-i xr 


0 «r 


0 0 


rv rc 


ro *-• 

cr c 

a < 


or it 


ro v 


CV vO 


r~> 90 


c tt 


IV tD 


M xj 


■W- 


CL tC 


- ■ in 




c in 


c . no 


H »| 


CM O 


cm m 


.0 rv 


<r 0 


M cj 


m «c 


rv 


O -r -1 


n •-» 


n CV 


^ CV 


-*-i CM 


n 


n CV 


n rr- 



O' 


CD 


















O 


CD 


















♦ 


♦ 


















Uj 


UJ 


















4 T 


9 


















O' 


O' 


















rv 


rv 
















♦ * 


r* 



















0 tr 

rS 

0 m 


‘0 ro 


0 ro 


0 0 


0 in 


O O 


0 .n 


0 rv 


0 rv. 

4 1 - 

c c 

• 

CD C 

• 

O CM 

* 

c rv 

• - 

0 O. 

• 

c. CC 

• 

' O tC 

• 

C M 

« 

c 

• 

C- c t 

— « w 

0 m 

to 

O CD 

lO 

0 0 

0 

0 0 

O 

0 Ot 

0 

0 0 

0 

O O 

0 

0 0 

0 

CD O 

0 

O 00 

Z? O 

(C *« 


e *-1 


a <*f 


c 


O 4 H 


c rl 


c rl 


C x- 


c: t-< 


C V* 

z z 

CM CM 


CM CM 


CM CM 


CM CM 


CM CM 


CM CM 


CM CM 


CM CM 


CM CM 


CM CM 


* * 


* * 


* * 


• * 


* * 


* * 


* * 


* * 


' • 


* * 



t) 


It 


it 


a 


11 


M 


u 


0 


it 




(X 


X 


X 


« 


ae 


X 


X 


X- 


X 




UJ 


U' 


u> 


Uj 


U’ 


Uj 


,u> 


u 


tv 




r~ 


♦— 


1 - 


f- 


►- 


b~ 


r- 


»— 


h- 




U: 


Ul 


u 


u 


U! 


UJ 


u: 


US 


u 




X 


X 


X 


X 


£ 


X 


X 


T. 


X 


* ♦ 


Mt 


9 


9 


9 


<4 


9 


-4 


9 


■c 



0 m 

QC 

O 

rr 

O O 

tr 

O XT 

cc 

0 O 

cr 

0 m 

X 

cd m 

X 

O -r* 

X 

0 ro 

X 

0 rv 

9 . *~ 

c • 

<t 

CD * 


c • 

9 

O • 

9 

O • 

•« 

c * 

>4 

c • 

4 

c • 

< 

c • - 

<r 

c • 

w >-» 

CM <0 

Ol 

O' to 

Q. 

M> 0 

a. 

<0 ro 

CL 

O CM 

OL 

CM O 

OL 

9 in 

Ol 

to tO 

X 

to rv 

X 

0 m 

a- a 

«L O' 


a: c 


aj « 


CO M- 


a 0 


O' -C 


O M. 


0 c 


0 c 


cr 0 


rv in 

O 

co 

O 

rv »v 

0 

rv m 

O 

rv cm 

O 

rv rv 

O 

rv n» 

O 

rv rv 

0 

TV 

0 

ID V 


m- 


c\ 

1 

ro 

• 

M 

1 

in 

r 

in 

t 

tC 

c 

< 

1 

rv 

1 

fv 



a . 


X 


a 


a 


ac 


a 


X 


X 


X 




4 

X 


o 

c 

Q 





o 


o 

o 



4 in 


O tv 


o 


to 


tf\ x> 


»v «o 



25 





4 . 2 NOLIN Sample Problems 

4.2.1 The sample problems for the NOLIN program have been chosen 
to exercise several options of the program. Graphite/Epoxy and 
Boron/Aluminum laminates have been modeled, employing material 
properties derived from the Air Force Composites, Design Guide. 

The initial problem is a +30° Boron/Aluminum laminate under 
uniaxial tensile loading. Ramberg-Osgood parameters were input, 
with exponents set equal to 3.0, and a maximum strain failure criteria 
was employed. The input data and the program output follow in 
sections 4.2.2 and 4.2.3 respectively. A graph of the axial stress- 
strain curve for the laminate is given in Figure 2 . Note that 
the laminate exhibits a non-linear stress-strain behavior from 
loading onset to failure. 

The second sample nonlinear laminate analysis problem is a 
0, +45 Graphite/Epoxy laminate under combined tension and shear 
loading. In this case the Ramberg-Osgood parameters were determined 
by the curve-fit routine from uniaxial stress-strain data supplied 
as input. The, input data and program output listing are given 
in sections 4.2.2 and 4.2.3 respectively. Output of. the. program 
is plotted in the form of axial stress versus axial strain for 
this case in Figure 3. 
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4.2.2 NOLIN SAMPLE PROBLEM INPUT CARDS 


C nOLIN UATA 


* 
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♦ 

★ 

+ 


NOLIN SAMPL.F PwnHLEf-'S 


VEHS ? MOD 3 


li 

2 . 


o, +45,-45 (.amIMATE:. FULL R-0 PapAM nUPVS*FlT, f/U- CHIT, 

♦ 30, -.TO LAmtmATE, R-0 PARAM INPUTT STRAIN Fail chit. 


spa t a 

NLAV=3. 

gll<l)s3*l .R02E+nS, 

£?2 f 1 > = 3 * 3 , 9 ont + o 4 i 

Vi2<i)a3*(3i?55&4on;- 

Gl2<l>=3*5.77f]!E«.0T; 

THICK<1>P0i 25»Q.5*0.25« 
TWETA(l) = A5>,,0,n,-45,» 
10PT«2 f EOPT«i.COpT«i, 


SlGllCl,l)? 6 , 89 F*fln.l 3 t 78 e *0 0 . 2 n. 67 E* 00 . 27 ,S 6 E*nO f i‘ 1 .455 + 110 , 41 1 3^6 0 ° 

48, 23 E-Q 0 , 55 . 126 + 00 , 62 . 016 * 00 . *8. 906 - 0 II. 

EPSl.l(l.l) s T,lt-o 3 -,?,^ 6 - 03 . 4 .E- 03 ^ 5 , 7 E- 03 . 7 . 6 P-n 3 ,l,E- 02 ,l, 31 E- 02 , 
SlG22<l,l5=6%?F4fln'!l3,78F*«0.?(I.67E + no,27,S6E*'iO,34,45F*OU,41,34F*00 

48,23t40n,S5,i?£4-on,' , '2.0TE*on, <sa. 90 £ + 0 n, 

1 ,65b-0P, 9 Ti 5E«P2. 2.8 4 F*0?. t 

SlG12tl,n s 6.89F4nn,l3 ? 78E*no,?n.6?F-0 0.27.S6F*^0.3 4,45F^00.41,3AF*On 



NOLIN SAMPLE INPUT CARDS CONT. 


46 1 23t-t-Ofi,'55..i?t- + on, # >2.ntt + on I 6r>.9nE't-OP, 

EPS12 < 1. 1 > s ? i 2 E,q 3 •. 5,26-0 3,3.6- 0 3'. 11. ,4|=-o 3, 1.5,2 = * 0 3,2. 6 - 0 ?.?, *2E*n2, 
3,30E“02»4,30F*-n2,«5.ftflEr'02i 

SOI 1 = *5 . QF*00, Sn??=-? ,5F«-0O» SOI 2 = 0 , 0 , 

IFC'Vs4, IPRINTsI, 

MATVPE(i)=l,l,i. 

Si*iT ( l)ai,32E03,s??T(l)B7,iftE + nt,Sl?(1 > = l7 O^E* Q2, FP UT< 1 )=6,63E-U3, 
EP22T<l) s 3 I 77E'!n3,f!AMA(l) = i 1 827Fs-n2,S1.ir(l) = 2 < 43E + 0' 1 !.'>22C<i) = ?,7t>F*O2, 
EPilC<l)»l«227fc-o?;Fp22r(i >31,456,02. 

STlFFafl ,100. 

a 12 ( 1 )s 3*^2.36?3Fin6, 

K$GM=50. SmL Ts 3, ft, ?T=lO,TNMT=2f 

$DATA 

NLAV=2, 

Ell<l> s 2*2»20E + 05, 

E22<1>=2M .246+05, 

Vj2<l>=2*O.(UE + 0Oi 
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THinK<l>=?*0,5Q, 

THET4 ( i ) = 30 . 0 n , 

lOPTal, 

STY<l>=2*1109,E*0O; 

SCY( i> = ?*i467 , fc*or>; 

TY(l) = 2 + 93. OE + 00, 

XM=3.O0,XN=3,0U. 

S011a50,» 
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SOl?s0 , , ' 

I F C N s 2 » 

SllT (1 ) = lj 00 ,0.S22T(1 >=103. ^ S. 1.2(1 ) =93 0 , FR1 1 T { 1 > =Q . 76" 02 , 
EP2?T(l> s ?»0E-0 2, Gama (15=3, 0E-.02 , S11C ( 1 ) =1.480 . Q , S22C ( 1> =160 , , 

EPj lC< l > S P t 7t-(j p, Ep-??C ( i. ) = o', 02. • 

STIFF=0,10U. 

KSGM=48 ,SkLT=4$ 


30 



NOUN SAMPLE 


O —* 
o -i a 

T -J U 


31 



data Input points for curve fit 
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lamina FAILURE CRITERIA 
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AT FIRST POST-FAILUPF LOAD POINT 
EXTERNAL APPLIEH 5 TRESS 
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♦♦♦♦ LAMINATE analysis interpolated to failure point ♦♦♦♦ 
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LAMINATK has failed* OUaO&ATtC INTERACTION FAtUJHk 

Quadratic = i* 03?0 for layer 2 

Quadratic a .9714 for l a y t r ? of previous load 
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APPLIED STRESS ANALYSIS 
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5» Computer Program Listings 


Source listings of both the UNI and NOLIN computer programs 
follow. The UNI program requires 20 K of computer core storage 
in a CDC 6600 machine, while the NOLIN program requires 60 K of 
core storage. No peripheral devices are required for either 
program for intermediate data storage. 
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ROUTINE ANGLE REDUCES ANGLES TO VALUES BETWEEN 0 AND P 1/4 FOR 
COMPUTING SIN AND COS 
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ROUTINE RESET SETS RELATIVEY SMALL STRESS TERMS EQUAL TO ZERO* 
TO AVOID CONVERGENCE DIFFICULTIES U.E.* VALUES ) HVAL) 
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ROUTINE rONVR CHECKS FOR CONVERGED SOLUTION DURING NfWTON-RAPHSON 
ANALYSIS, ALSO CHECKS ITERATION LIMIT ANO DIVERGENCE LIMIT, 


+ 

ft- 

\ 

ro 

tv 


ft- 

o 


o 

*> 


U- 
























•-» 
























O K“ 
























- -i 
























- i 
























~ 00 






• 
























co 


















o k 






3 c 
























— no 


















— ■ 






»-* ro 





in 













k- 






IT 





ft- 


© 











C 






in c 



© 


© co 













W Z 






— j_ 

•> 


* 


in 


ro 




o 







*■ 






to 

no 


ro 


ro o 






O 







- o 






X o 

3 




i- 


c 




M 







— X 






<t ID 

— 


o 


o 


h- 




• 







» a 







r-i 


K 


ft- O 













C D 






• « 

ID 




ID 


O 




X 







\C * 






~ © 

in 


o 


O 


ID 




. X 







~ lA 






~ UJ 

■N. 


ID 


ID 






ro 

-i‘ 






id a 






r-t O 

CC 










« 

3 






m u; 






» • 

3 


«-* 


O 


C 




(V 

C 




Q 







u 

ro o 

in 


in 


— X 


X- 




X 



in 


X 


\ 





UJ 

3 • 

— 


a 

* 

o a 

u 

a 




rr 



a 


.a 







X 

m — o 

in 


UJ 

o 

• 3 

UJ 





» 

U’ 


U! 


3 


*-* 





o 

1- O U) 

(X 

X 

• 

UJ 

UJ • 

I 

• 




© 

o 






»- 






_i in • 

< 

3 

i- 

X 

Z UJ 

o 

UJ 




© 

z 






UJ 



co 



» 


in 

O 

o 

• 3 


-i 





UJ 


© rvj 


© OJ 


in 





UJ 

^ in ~ 



• 


• 


• 





ID 


oj rv 


ro rvi 


V 



> 


a 

n cc • 





b- — 

UJ 





© 

X 


r- r- 


r- ft- 


z 

r— 


< 


2; 

no « ro 

0 IT 

u 

IT 

z 

►- rn ft, «j ir 

u 

r. 

if 

rv 


c 

u 


^ — u 


r— m* 


o 


> 

* 


UJ 

“5 -> 

ro 

Ui 

U. 1 O ft 

o 

U‘ 1 3 r- 

z 

3 

ft- 

u 

UJ 

If 

> 

u 

* • a 

u 

m- m 

u 

*— * 


<* 



ID 

— • 

~C0 3 - 

3 — ■ <o 

►-3 k- ro 

UJ 


ro 

3 

3 

■— 

z 

3 

x x> r' 

3 

sC *> 


in z 

n 

-J 


r*4 

X IT ~ 

fO 

Z ro 

z u 

ft- 

Z ~ U U u 

<D 

u. n 


Z II 

Z 


o 

z 

— ■ »-* 

Z 

«—> 

z 

z o 



n 

a 

UJMI 0 

30 

i- 3 

*-♦ »-» o 

<1 

*-* z *-> c 

cr 

►-« 

o 

»— i 

*— 

o 

u 


UJ Ui o 


UJ UJ 


UJ X 

z 



z 

> 

no in 

— b~ 

K — ’ 

ft- C H 

a 

ft-'-'ZCZft- 

UJ 

c z 

K 

K 2 

K 

ft- 

1 

K 

ft- b~ b- 

1- 

ft- ►- 

1- 

X 1 

c 

II 

ro 

ft- 

Z 

cc* w 

U 

Z U- 

z - 

u* 

z - o - c 

> 

— o 


Z O 

Z 


z 

z 

i— < •— * 

z 

ft- 1 *-« 

z 

*-* o 

o 


>- a 

o 

O 3 U 

•-* o o *-« O U- o 

i-ou.ou.oo 


u u 

o 

COO 

o 

o 

o a a: c 

o a a 

o 

c u 

Hi 

z 

-i 


o c in»— 

OID O C 

o e 



o 

M #— » 

c 

o •- 

O 

ir 

z 

o 

» i ir 

c 

J X 

o 









o 

in 


© 


© 


© 

IT 



CM 




IT 









ro ■ 

ro 




in 



k- 



CC 


X 


O' 









ro 

ro 


*o 


ro 


ro 

ro 



ro 


ro 


ro 


o 




O o 





o o 


u o 


o 

a 


o 

oo 







in 

a 

LU 

> 

Z 

3 

T 


o o o «— « cm ro 
«— < 

o o o © © © 

e e c c c c 

o © o o o o 

D o © O © © 


4 O’* \C^W(V ^ <0 O O (VI ft « 

~ IV ruro < <t 

c c c ccooooo oooocc 
c c c c c o c c o c o c c c e c 

C CO ooooooo © © © © © © 

©©©o©ooo©o © © © © © o 


o co Mf 
IT IT IT 

o o o 
cc c 
oo c 
© © © 


U' IT ^ 
If' If If 

o © © 
o © © 
o o o 
0 0 3 


© 


r-f-ft-c.c©- 4 f\j 

OO 

eccccocc 
cooc coco 
oooooooo 


88 



r* 

\ 

n 

cm 

N. 


O 


x 

> 

z 

o 

o 


• 

I 

m 

sO 

m 


UJ 

> 

uj 

_J 

x 

in 

a 

• 

i 

(V 

z 

o 

m 

a 

UJ 

> 

z 

3 

X 





a 



0 


ro ^ 






-* 


»-» • 

• 


• Z 



UJ 


• » 

>- 


• ' 



l D 


♦ >* 



— rvj 



a 


X 

_J 


—« ♦ 



UJ 


z 

UJ 


© t-> 



> 


O X 

X 


in — 



z 


*-• CD 

• 


• u. 



o 


>- CO 

* 


*» M 



o 


< * 

X 


^ Q 




a 

X * 

CD 


» m. 




* 

UJ X 

» 


z — 



J- 

* 

>- a> 

* 


* z 



c 

in 

M » 

X 


fVI ♦ 



z 

UJ 

— * 



♦ !-• 




- ID 

• >- 

-1 






m x 

• 

UJ 


- u. 



in 

• UJ X x 

X 


CD *-• 



u> in > 

GO CD 

• 


in o 



o 


cm in 

m 


» » 



— o 

Ui O 

•> 0 

X 





in 

m 

♦ * 






• 

• 

^ X 






m in 

• in 

0 GO 

» 


z u. 



— in 

►- in 

• » 

o 


4 



UJ UJ 

CD UJ 

ro « 

>- 


►- o 



co a 

• a 

*-• X 

X 


» 



•> t- 

»- 

» 



© — 



~ in 

in 

♦ X 

z 


in — 



X 

a 

CO 

o 


» » 



-<* 

o 

z in 

in 


- z 



• a 

a x 

o a 

* 

a 

«-> * 



in o 

x a 

*■ 

•> 

i- 

» CM 



• u 

UJ U- 

l- X 

X 

•-* 

p-» ♦ 



m 


<3 -(J- 

CO 

z 

•— • 





X * 

• 

m 

O - 



UJ z 

UJ z 

UJ « 

* — 

*— 

in — i 



ro o 

> o i- x 

V \ 

*—* 

• o 



*— < 

»— i #-• 

»-» UJ 

* 

— 2 

*-« in 



cm 

*- ►- 

— > 

z >- 


>~ 



• 3 

< 3 

o < 

CD X 

• — — 

<i ^ 



x _i _j _l 

• _J 

in 

JO 



— O 

UJ o 

X 

* -i 


♦ n 



* in 

a m 

00 0 

• UJ 

»-» F* f*~ 

— in 






X X 

2 —i — 

H — 



* 

* « 

X 

a o 

♦ • 

* U* 




^ — * 

w w 


>0 *> 

*> 3 







a 

r- w z m 

z 

CM Z 

»- h~ 

i- >- 

»“ 1- 


uj u: cr uj •- a 

a 

a a 

<z «x 

<c <C 

< <r 


a k k 

r> k- z 

3 

2 3 

2 z 

z x 

X X 


K H H 

»-« 2 K 

H- 

►- *- 

a a 

a a 

XX 


*-• a cr 

ox o a 

uj a UJ 

o o 

o o 

o o 


2 3 3 

C 3 O 1- 

a 

p-i a 

u u 

u u. 

u u 









^ CM 


r- 


© © 

© © 

CM © 

-4 CM 



O' 


o o 

in cm 

cm m 

<* <f 



ro 


IT 

in i" 


r- 







•— » ^ 

^ F - 



o u 


C 

z 

UJ 


ro co r- x* 

— r- 

©O—I 

CM CM 

CM CM 

CM CM 

CM 

cv cv co n <t 

4 4 

in if; in 

in if i 

if. in 

If. in 

M. 

t~* —• •«- ^4 

CM CM 

CM CM CM 

CM CM 

CM CM 

CM CM 

CM 

c c o c c 

© © 

CO© 

© © 

© © 

© © 

e 

© © © © © 

© © 

© o © 

o © 

© © 

© © 

© 

© © © © © 

© © 

© © © 

© © 

© © 

© © 

© 


89 



ROUTINE LAMTST PERFORMS FAILURE ANALYSIS* STIFFNESS TEST » ULTIMATE 
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TOMPUTFS quadratic failure criteria coefficients 
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FUNCTION S'JBHROGHAM SINT dfrtimines interpolation factor for 
STPPS OR STRAIN FAILURE 
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PUOCK OATA SUBPROGRAM INITIALIZES VARIABLES 
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